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Abstract. This paper presents a method to find new De Bruijn cycles based on 
ones of lesser order. This is done by mapping a De Bruijn cycle to several vertex 
disjoint cycles in a De Bruijn digraph of higher order and connecting these cycles 
into one full cycle. We characterize homomorphisms between De Bruijn digraphs 
of different orders that allow this construction. These maps generalize the well- 
known D-morphism of Lempel [13j between De Bruijn digraphs of consecutive orders. 
Also, an efficient recursive algorithm that yields an exponential number of nonbinary 
Dc Bruijn cycles is implemented. 



1. Introduction 

The Lempel homomorphism between binary de Bruijn graphs of consecutive orders 
has been used by many authors for over thirty years to construct new de Bruijn se- 
quences using a given one of smaller order as in [3], [5], [13], or to obtain results about 
the linear complexity of binary sequences as, e.g., in [1], [8j and references therein. 

Ronse gives an attempt in [16] to generalize Lempel's construction by presenting a 
single function between two non-binary De Bruijn digraphs of consecutive orders that 
has a similar effect. It appears that there is hardly any other attempt in literature to 
generalize this construction despite the importance and simplicity of the approach. It 
is worth mentioning that [H] presents a collection of homomorphisms from higher order 
De Bruijn graphs to lower order ones, but these do not nearly enjoy the properties of 
the Lempel homomorphism as will be shown in the sequel, and thus can not be used 
to construct De Bruijn cycles. In this paper we generalize Lempel's homomorphism 
by describing and characterizing a class of homomorphisms between two De Bruijn 
digraphs of arbitrarily different orders but with the same alphabet, the direction of 
these functions being of course from the higher order digraph to the lower order one. 
That is, we both generalize to non-binary alphabets and consider De Bruijn digraphs 
of non-consecutive orders. The definition of the proposed class is exceedingly simple. 
A homomorphism belongs to it if any cycle in the lower order digraph is the image of 
a fixed number of nonintersecting paths in the higher order digraph and of the same 
size as the given cycle. This fixed number depends only on the difference between the 
orders of the digraphs. 
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Besides the fact that these homomorphisms attest to the structural similarity among 
different order De Bruijn digraphs, their significance stems from their applicability to 
interesting problems. For example, noticing that a generalized feedback shift register 
sequence (GFSR), described for example in [15], is indeed a cycle in an appropriate 
De Bruijn digraph of some large alphabet, the recent work P exploits the proposed 
homomorphisms to obtain a virtually endless number of parallel streams of pseudo- 
random numbers based on a given GFSR. The streams are simply the inverse paths of 
the GFSR sequence. These streams can be made very much "uncorrelated" by selecting 
a homomorphism that is highly nonlinear and two De Bruijn digraphs with a large 
difference between their orders. Furthermore, this is done efficiently by eliminating the 
traditional preprocessing effort needed for jumping through a fixed stream of pseudo- 
random numbers. 

Another important and relevant problem is the decomposition of a De-Bruijn graph 
or digraph into a set of disjoint cycles called a factor, see \X0\, [13j for instance. By 
definition of the proposed homomorphisms, it can be seen that a factor in the lower 
order digraph readily induces a factor in the higher order digraph. 

In the current work we only deal with the particular case when the vertices of the 
lower order digraph are all on one cycle, i.e. on a De Bruijn cycle, the inverse image 
paths turn out to form a number of non-intersecting cycles that include all the vertices 
of the higher order digraph. Thus cross-joining these cycles into one big cycle produces 
a De Bruijn cycle in the higher order digraph. The focus of this paper is to describe and 
characterize the proposed class of homomorphisms and to present particular members 
that simplify the task of cross-joining those inverse images into a new full cycle. 

In Section[2]we include definitions and known results that are useful for the rest of the 
paper. In Section [3l we present an immediate generalization of Lempel's construction 
to the case of nonbinary alphabet size q. This case turns out to be much richer than 
the binary it allows more homomorphisms and more alternating strings to be 

used for the construction. The number of resulting De Bruijn cycles of order n is 
exponential with a huge base of size q(f)'^{q) where (p is the Euler totient function. We 
describe two variants of an algorithm that performs this cross-join construction. We 
also present an efficient implementation of this algorithm that is based on a formula 
that calculates the location of the cross-join. 

In Section H] we introduce the class of homomorphisms discussed above. Propo- 
sition 14.11 characterizes all De Bruijn homomorphisms while Theorem 14.31 provides a 
characterization of our proposed class. Proposition 14.51 studies a special homomor- 
phism between two binary De Bruijn digraphs with an order difference of two, and 
whose application yields two cycles of different lengths that can be joined relatively 
easily into one full cycle. 

2. Preliminaries and Basic Results 

2.1. Terminology. For positive integers n and q greater than one let be the set 
of all g" vectors of length n with entries in the group Zg of residues modulo q. When 
the group structure is not needed we will sometimes refer to elements of this group as 
symbols. We will interchangeably use the terms vector, string and word in the sequel 
to mean a sequence of consecutive symbols. An order n De Bruijn digraph, Bn{q), is 
a directed graph with as its vertex set and for two vectors x = (xi, . . . ,a;„) and 
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y = ...,?/„), (x; y) is an edge if and only if yi = x^+i; z = 1, . . . , n- 1. We then say 
X is a predecessor of y and y is a successor of x. Evidently, every vertex has exactly 
q successors and q predecessors. Furthermore, two vertices are conjugate if they have 
the same successors. 

A cycle in Bn{q) is a path that starts and ends at the same vertex. It is called 
vertex disjoint if it does not cross itself. Two cycles or two paths in the digraph are 
vertex disjoint if they do not have a common vertex. The weight W{w) of a word 
w is the sum-in Zg-of all elements in w. A translate of a word (or cycle) c is the 
word (or cycle) c + A where A is any scalar and addition is component-wise. A cycle is 
primitive in -B„(g) if it does not simultaneously contain a word (xi, . . . , x„) and any of 
its translates. For a positive integer r a function d : — * is said to be translation 
invariant if d{w + A) = d{w) for any w G Z^ and scalar A. 

A De Bruijn cycle of order n is a Hamiltonian cycle in Bn{q), i.e. a cycle that 
visits each vertex exactly once. Alternatively, a "linear" De Bruijn sequence of order 
n is a sequence of symbols from Zg in which every vector of Z^ occurs as a substring 
exactly once. For example, 0001110100 is a binary De Bruijn sequence of order 3. Its 
corresponding De Bruijn cycle is 000 ^ 001 Oil 111 110 101 -> 010 
100 — > 000. For the rest of this paper we will denote such a De Bruijn cycle-and in 
fact any cycle-as [11101000], and consider any rotation of it to be an equivalent cycle. 
We will denote the linear sequence as (0001110100) and will simply write i?„ to denote 
the binary De Bruijn digraph of order n. Finally for the significance and many known 
algebraic, combinatorial and graph-theoretical methods of construction of De Bruijn 
cycles we refer the reader to [9], |T0], and references therein. 

2.2. Lempel's Homomorphism. Define a map D : B„ B„_i by 

D{ai, . . . , a„) = (fli + a2, 02 + • • • , a„_i + a„) 

where addition is modulo 2. This function defines a graph homomorphism (see below) 
and it is known as Lempel's D-morphism due to the fact that it was studied in |13j . 
although it can be traced back to [12]. Note that -D(x) = D(x + 1) for all x G Z2 . We 
define the dual of a cycle c, to be its bitwise complement c -|- 1. c is called self-dual if 
it is a rotation of c + 1. The following facts are proved in [T3] . 

Fact 2.1. A cycle of length p in Bn-i is the D-morphic image of two primitive, vertex 
disjoint cycles of length p in Bn if and only if it has an even number of ones. 

Fact 2.2. A cycle of length p in Bn-i is the D-morphic image of a self-dual cycle of 
length 2p in Bn if and only if it has an odd number of ones. 

Two cycles are called adjacent if a vertex on one cycle has a conjugate on the other 
cycle. Swapping the successors of these two conjugate words joins the two cycles into 
one larger cycle. This is why any pair of conjugate words is called a cross-join pair, (a 
similar concept of a cross-join tuple can be defined for q > 2.) By Fact 12. H if c^-i is 
a Hamiltonian cycle in B^-i, then every word in Z2 is either on c„ or on c„ + 1, the 
two primitive pre-images of c„_i by D. Lempel used this idea to construct De Bruijn 
cycles recursively by rejoining c„ and c„. The most obvious cross-join pair is the two 
alternating strings of size n 2;^ = 010 . . . and its complement Zn which can not be on 
the same cycle. Recently, this method was implemented in [3] with an efficient, linear 
code, and more recently done even more efficiently in [5] with a jump from a given De 
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Bruijn cycle in to a higher order -B„+fc, for some integer k that is a power of 2, by 
pre-computing the effect of an iterative apphcation of the D-morphism. 

3. Generalization To Nonbinary Alphabets 

The discussion in this paragraph aims at an immediate generahzation to the q-arj 
case of the recursive construction of full cycles using the Lempel homomorphism as in 
[3|,[5]. In fact, p£j gives an attempt to generalize this homomorphism by considering 
the single function defined by (p{xi, X2) = X2 — X1. Note that this function is translation 
invariant. 

3.1. Recursive Construction. We begib with the following proposition which relates 
translation invariant homomorphisms to primitive cycles. It generalizes Fact 12.11 and 
Fact 12.21 by finding all translation invariant functions that can be used to mimic the 
Lempel construction of de Bruijn cycles in the g-ary case, as described below. 

Proposition 3.1. Let T be a vertex disjoint cycle of length I in and -D„_i,i be 

the homomorphism from Bn{q) to defined as 

Dn-1,1 = idi{Xi, X2), di{x2, X3),..., di{Xn-l,Xn)) 

where di{xi, X2) = axi + (3x2, a and (3 are elements in Zg such that gcd(/3, q) = I and 
a + (3 = 0. Then 

(a) W{r) = if and only ifV is the image by -D„-i,i of a primitive cycle of length I 
in Bn{q). 

(b) W(T) = X ^ if and only ifV is the image by Dn-1,1 of a cycle 

(3.1) C = c ■ (c + f3'^X) ■ (c + 2f3'^X) . . . (c + (r - 1)(3-^X) 

obtained by concatenating c with its translates, where r = q/ gcd(A, q), c = [xi, . . . ,Xi] 
is an appropriate primitive cycle of length I and j3~^ is the inverse modulo q of (3. 

Proof. First, the fact that -Dn-1,1 is a homomorphism is straightforward. The reason 
for the index notation in -Dn-1,1 and di will be clear in the next section. We will 
now prove (b). Let F = [71,..., 7^] be an image of a cycle C in Bn{q) given by 
(13. ip . Then 71 + . . . + 7^ = {axi + (3x2) + . . . + {axi-i + I3xi) + {axi + (3xi+i) = 
a{xi + . . . + xi) + (3{x2 + . . . + xi+i) = (a + I3){xi + . . . + xi) + f3{xi+i - xi). Since 
a + j3 = and is the first element of (c + P~^X) it follows that the last expression 
is A. Conversely, let F = [71, . . . , 7/] be a cycle in Bn-i{q) with weight A 7^ 0, be an 
arbitrary value in Zg and define c by letting 

Xi = /?~-^(7i + . . . + 7j_i) + xi, i = 2, . . . , /. 

Evidently r is the smallest value such that c + r(3~^X = c, hence C as defined in (13.11) is 
an rZ-cycle in Bn+i{q) whose image by Dn-1,1 is F. To see this note that = —1 and 
for i = 1,. . .,1-1, axi + (3xi+i = a/5"^(7i + . . . + 7i_i) + + (71 + . . . + 7i)+/5xi = 7^. 
Also for j = 0, . . . , r - 1, axij+i + (3xij+i+i = a{xi + j(3~^X) + (3{xi + (j + l)(3~^X) = 
axi + ia[3~^X + (3xi + (j + 1)A = axi + (3xi + X = a/5~"^(7i + . . . + 7;_i) + aa;i + /9a;i + A = 
— (71 + . . . + 7i_i) + (a + l3)xi + A = 7/. The last equality follows from the definition 
of A and the assumption a + (3 = 0. Moreover, note that for j = 0, . . . ,r — 1, 

the noncychc sequence (x/j+i, . . . , x^j+i, . . . , x/(j+i)+„) does not contain any 
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translates, otherwise two translates would have the same image by -D„_i,i violating the 
assumption that F is vertex disjoint. 

Part (a) can be done similarly by defining r = 1 and = xi in the proof of 
(b). □ 

Technically, Part (a) is a special case of (b). We separate them for easy reference 
below and because Part (b) straightens Proposition 98(b) in [16]. The choice of a = 
q — 1 and (3 = 1 corresponds to Rouse's function. The algorithm below shows how the 
above proposition helps to construct q^p{q) De Bruijn cycles of order n with alphabet 
size q using one De Bruijn cycle of order n — 1 and a homomorphism -D„_i,i satisfying 
the requirements of Proposition 13.11 A definition and a lemma are in order. 

Definition 3.2. For a given value A G Zg we define an alternating sequence 9^^^ to 
be the infinite sequence {ei, 62, 63, . . .} such that ei = and Cj+i = Cj + A for i > 1. 
Let oi^^ be the suffix of size n of O^^'^K We also define the cycle An to be the self loop 
(A ... A; A ... A) in Bn{q), for any digit A. 

Note that the iterates {iX + Oi^^ —>■ {i + 1)A + ^n^"*); z = 0, . . . , g — 1 make a cycle of 
size q in Bn{q) for any A that is coprime with q. 

Lemma 3.3. For any element 7 7^ m the inverse image of the self loop 7n-i (in 
Bn-i{q)) by Dn~i,i of Proposition [XTl is the cycle of length q in Bn{q) formed by the 

alternating string (fn '^'^ ■ 

Proof. Let [0,X2, ■ ■ ■ ,Xn] be the inverse of 7n~i that starts with xi = 0. We need to 
show that Xi = {i — l)/?~^7 for i > 1. This is satisfied when i = 1. Proceeding by 
induction, assume this is true for i. Then aXi+Pxi^i = 7 implies {i — l)aP~^j+Pxi+i = 
7. Since aP~^ = —1 it follows that Xj+i = iP'^j- Finally, since both P and 7 are 

coprime to q so is the product /3~^7. Thus the cycle formed by 6^ "'^ actually has 
length q. □ 

The following algorithm is a direct generalization of Lempel's binary construction 
outlined in Subsection 12. 2[ 
Algorithm A 

Given: a and A G Z^, (A coprime with q), a homomorphism -Dn-1,1 with a and (3 
satisfying Proposition 13.11 and a De Bruijn cycle F of order n — 1 written linearly as 

( A, . , A , 7^, . . . 'Jqn-l , \, . , . 
n—1 n— 1 

1) Form the inverse image {xij}, j = 1, . . . , q^~^ + n of F that starts with a. 

2) For i = 2, . . . , g let {xjj} = {xi^j} + {i - 1)(3-^X. 

3) Regarding the sequences formed above as cycles in Bn{q), erase the edge that emerges 
from the first n-vertex in all cycles except the one that starts with a + (g — 1)/3^^X. 

4) Erase the last edge (that reaches back to the first string) in all cycles except the 
starting cycle {xij}. 

5) For i = 1, . . . ,q — 1, create an edge from the string (xi,i, . . . , Xj,„) to the one that 
starts with (xi+1,1, . . . , Xi+i,„). 

6) For i = 2, . . . ,q, create an edge from 

(Xi^gn-i , . . . , Xi^gn-i^„_l) to (Xi_1^27 2^4-1,3! ^i-l,n+l)- 
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There are q inverse images of F that collectively include each string of size n exactly 
once. The algorithm outputs one full cycle of order n obtained by connecting the former 
q cycles. Since -D„_i^i is translation invariant, exactly one translate of O'^ '^^ occurs on 
each separate cycle. By Lemma the first n- vertex in each inverse cycle is a translate 
of^r^'^ Be ginning with 9^ + a. Step (5) iterates through all the other translates 
oiOn the last one being + a + {q — l)(3~^\. Since (xj+i j^n-i, a;j+i gn-i+„_i) 

and (xj^i, . . . , Xi^n) are both predecessors of (xj+i^i, . . . , it follows that the vertex 

(xj_2, • • • , a^i^n+i) is a successor of (xj+i_gn-i, . . . , Xj+i^gn-i+„_i), this justifies step (6). 
Steps (3) and (4) of course do the necessary deletions before individual cycles are 
joined in steps (5) and (6). 

Figure 1 illustrates the above construction for g = 5, n = 2 and a = where each 
row stands for an inverse cycle and the first row represents the cycle in Step (3), the 
alternating string is viewed as the first vertex in each cycle and the string {x + i,i,i + l) 
as the last vertex. 
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(3,4,0) 




^(4,0,z+3) 


^ (x+4,4,0)' 




(4,0,1) 




(0,l,z+4) 



Figure 1. 

Proposition 13.11 (b) can also be used to construct the same q^{q) cycles that are the 
output of Algorithm A. Note that, regardless of the value of g, F has weight zero. Let 
F~ be the same as F with exactly one single occurrence of a letter 7 deleted from the 
only run of n — 1 7s, where 7 is chosen such that gcd(7, g) = 1. Then the length of 
F" as a cycle in i?„_i(g) is g""^ — 1 and gcd{W{V~),q) = 1 and thus constructing 
the inverse image by -Dn-1,1 of F~ yields a cycle in Bn{q) of length g" — g. Using 

Lemma [3.31 we see that the inverse image of [7] is the cycle formed by iterating 6^ ^\ 
This shows that the g vertices of Bn{q) that are not in C~ form exactly one cycle 
which is, once cross-joined to C~, forms a full cycle. This construction is summarized 
as follows. 

Algorithm B 

Given: a and A G Zg, (A coprime with g), a homomorphism -D„_i^i with a and (3 satisfy- 
ing Proposition l3.1l and a De Bruijn cycle F of order 1 written as [ A, . .j , A^ , 7„, . . . 7g"-i] 

n-l 

1) Let F^ be [7„, . . . , 7gn-i, A, • - ; -^J and let F = (71, . . . , %n-g) be a concatenation of 

n-2 

g copies of F~. 

2) Let Xi = a. 

3) For i = 2, . . . , n + g — 1 let Xj = + P'^X. 
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4) For i = n + q, . . . ,q'^ choose Xi such that Q;Xj_i + = ^i-n-q+i- 

Finally, to see that the outcome of this algorithm is the same as that of Algorithm A 
the key idea is to note that, for a given a and A, the outcomes of both algorithms start 
with the q alternating strings in sequence followed by q sections of length g""^ — 1 each 
that are all translates of each other. 

3.2. EfRcient Implementation. The two algorithms above are mainly theoretical, 
establishing the possibility of the construction. To actually implement a recursive 
construction efficiently, note that we need to keep track of the locations of the constant 
strings within the De Bruijn sequence being inverted, since the inverse of a constant 
string is an alternating string where the cross-join must take place. A De Bruijn 
sequence that results from input i and A in Algorithm A will be said to be of type 
(i; A). Following the notation of Annextein [3], we say that a sequence representation 
of a cycle is oriented at 0„ if the all zero pattern is written at the end of the sequence, 
with the understanding that 0„+i is the first vertex in the cycle, e.g., [11101000]. With 
this representation of a cycle C, |x|c denotes the index within C of the ending symbol 
of x, provided that x is a substring of C . 

Proposition 3.4. Given a homomorphism D, let r„ he any order-n De Bruijn sequence 
oriented at 0„, where n>l. Let Cf, i = 0, . . . ,q — 1 denote the cycle started at symbol 
i that is the pre-image ofTn by D. Let tI^'_^1 be the order-(n+l) De Bruijn sequence 
obtained via Type {i; A) . Then for 7 G Z^, 7 7^ 0, 

{(g — m)g" + m i = 

(m' — m){q^ — 1) i 7^ and m < m' 

{m' — m){q^ — 1) + i 7^ and m > m' 

where m' = (— z)A~^ and m = (7 — i)X''^, all calculations being in Z^. 

In the proof below, phrases like higher cycle, next cycle, and top cycle refer to the 
diagrams in Figure 1 and Table 1. 

Proof. Let the homomorphism be induced by the function d{xi,X2) = oix\ + (5x2- The 
constant n-string /?A, denoted (/?A)„, has an inverse image in each cycle Q. We denote 
by *6'^+i the one that belongs to C,. To prove (i) note that there is an element m such 
that 7 = mA for some m, since A is relatively prime to q. Thus we have m = 7A~^ and 
7„+i resides on the cycle Cm\- Furthermore, it is the last string of size n + 1 on this 
cycle because Co is oriented at 0„+i. 

Now, Type (0, A) uses the first |°6'£^^|co and then proceeds through the alternating 
strings by adding one digit from each of the cycles Ca, C2A, C^a ("^ digits) and 
down to the last cycle. 

All the (g — m — 1) cycles that are visited after Cm\ is first traversed are used in full, 
thus adding {q — m — l)g"' digits to the length of |7n+i| {o;a). 

^ n+1 

Finally, the second and final visit to Cmx happens just after the location of the alter- 
nating string ei^l^. Since |°6'J,+i|c„;, = \^Onli\co, we use g" - \°6i+i\co digits of Cmx 
after this visit. By adding all the quantities together, we get 

|7n+i| {o;A) = |X+ilco + m + {q-m - l)g" + g" - |°6'J,+i|co = m + (g - m)g". 

n+1 
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Table 1. (a) Type (0;A). (b) Type (i; A) with /a(^,0) < lx{i,-f). (c) 
Type (z; A) with /a(^,0) > 



(ii) When i 7^ 0, type {i; A) is obtained by starting with the cycle Cj, as in Table 1 (b) 
and (c), and connecting through the alternating strings O^li- By definition, |7n+i|p{o;A) 

^ n+l 

is the distance from the start in a cycle oriented at 0„+i. For j G Zg, let the number 
of cycles that separate Cj from Ci through steps of size A be denoted lx{i,j). It is 
easy to see that lx{i,j) = (j — We will consider two cases. First, assume 

/a(^,0) < lx{i,'y). Then since Cj is the top cycle, the start of Cq is encountered during 
the upward traversal. Starting the count at the beginning of Co, we scan |°^i+i|co ~ 1 
digits, i.e., just before the end of the alternating string. The next set of digits are 
from C_A, starting just after the end of the alternating string till the end of this cycle, 
thus using (g" — |~'^^i+i|c_A) digits. To reach the constant string 7„+i we need to 
repeat the previous step for all /a(^,0) cycles above Cq. Once this is done, the start 
of Ci is reached. So then we reach the end of *6'^+i, with |*6'^+i|ci digits and go down 
along the last digit of each alternating string, down to cycle Cj+(g_i)A, this adds (g — 1) 
digits. Next, •jn+i is located at the end of cycle C^ so we need to scan upwards the 
(g — Ixihl) — 1) cycles underneath C^, using (g" — 1) digits of each. Finally, we scan 

the last (g" — P6'^+i|ct,) digits of C^ to reach 7n+i- 
Hence, it follows that |7„+iL(i;A) is the sum 

^ n + l 

(/a(^, 0)(g" - 1)) + {\¥^'Uc?l +{q-l) + {{q - h{^, 7) - - 1)) + - K%\c,) 
Using the fact that p^|,^i|c = |(/3A)„|r„ for all j, this simplifies to 

(ZA(^,o)-/A(^,7))(g"-i) + g"+'- 

The second case is Ixih^) > h{hl)- Then C^ is above Cq, see Table 1 (c). Starting 
at the beginning of Cq we reach the end of C^ by scanning upwards [q"^ — 1) digits of 
the (/a(^,0) — lx{i,'y)) cycles. This completes the proof. □ 

The most important feature of Proposition [33] is that the distance |7„+i| (i;A) depends 

-'- n+l 

only on the parameters 7, A and i but not on any previously constructed cycles. The 
homomorphism parameter f3 is only needed to determine which constant string in F^+i 
is the image of ^fj+i- 

Now we can write down a pseudo-code that is capable of constructing an exponential 
number of order n De Bruijn cycles for any alphabet size g, based on the trivial low order 
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cycle. We will use the order-1 cycle Fi = [12 . . . (g — 1)0] for odd q. Note that the weight 
l^(ri) = if and only if q is odd, thus allowing the construction in algorithm A. Note 
that the intermediate cycles can be 'lifted' to higher orders via distinct homomorphisms. 
Hence the parameter (3 can be changed in 0(g) ways. Since there are g0(g) type (i; A) 
cycles for a given homomorphism, we can construct (g0(g)'^)"^^ cycles of order n using 
Fi. The following algorithm is an implementation of Algorithm A above. It takes as 
input an alphabet size q and three g-ary strings B = {132, ■ ■ ■ , Pn), L — (A2, .... A„), 
/ = (^2, • • • , ?'n)- For each recursive step j, j3j determines the homomorphism to be 
used for this step, while (i^; \j) determines the type of cross-join. 
Algorithm AA 

Input: g, and three strings B = {(32, ■■■,f3n), (A2, . . . , A„), 7 = (^2, . . . , i„) 
Output: a unique order-n De Bruijn cycle F„ 

(1) Let zi = Ai = 1 and fJ^'^^ = 12 ... (g - 1)0 

(2) For j = 2, . . . ,n repeat the following steps: 

(3) Let a = q — Pj and calculate Co as the inverse of r^^-i'^^-i^ by the homomorphism 
induced by d{xi, X2) — olxi + l3jX2 

(4) Construct Ck = Cq + k;k — 1, . . . ,q — 1, [component- wise addition by k] 

(5) Let 'y = Pj X \j{ mod g) 

[Steps (6)-(10) determine the location pos of the cross-join within each cycle Ck] 

(6) If j = 2 then let pos = 7 and skip to Step (11) 

(7) Let m = (7 — ij_i)Xj\{ mod g) 

(8) If ij-i = let pos = (g — m)q^~'^ -\- m and skip to Step (11) 

(9) Let m' = {—ij-i)Xjli{ mod g) 

(10) If m < m' let pos = {m' — m){q^~'^ — 1) 
otherwise pos = {m' — m){q^~'^ ~ 1) + 

[Steps (11) to (23) connect Co, . . . , Cq_i by the string 9^] 

(11) Let newstring be an empty string 

(12) If ij = do steps (13) to (16) otherwise do steps (17) to (22) 

(13) Append the prefix of size pos to newstring 

(14) For A; = 1 to (g — 2) append the pos*'* digit of the cycle Ckx^ to newstring 

(15) Append the last (g-'"^— pos-|-l) digits of C(5_i)a^ to newstring 

(16) For k = {q — 1) down to 1 

Append the first (pos -1) digits of Ckx^ to newstring 
Append the last (g-'^^— pos) digits of C(^k-i)\j to newstring 

(17) Let m' — —ij\J^{ mod g) 

(18) For A; = to m' - 1 

Append the first (pos -1) digits of C^kXj to newstring 
Append the last (g-'"^— pos) digits of C(^k-i)\j to newstring 

(19) Append the first pos digits of C,;^, to newstring 

(20) For A; = 1 to (g — 2) append the pos*'^ digit of the cycle Ci^+kXj to newstring 

(21) Append the last (g-'"^— pos-|-l) digits of Ci.^(^q_i)\. to newstring 

(22) For A; = (g — 1) down to (g — 1 — m') 

Append the first (pos -1) digits of CkXj to newstring 
Append the last (g-*"^— pos) digits of C(^k-i)\j to newstring 

(23) rf''^'^ = newstring. 
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An implementation of this algorithm in C++ is available online, see |2]. We chose 
to use odd q because the position, pos, of the cross-join becomes especially simple to 
evaluate, see step (6). For even q, only the base case Fi has to be changed to some 
order-two cycle F2 and the position of the constant string 72 has to be located by 
searching the string F2. Table 2 shows all possible output that can be constructed 
from Fi = [120] for g = 3. 









2 


(1,0,1) 


[120221100] 


(2,0,1) 


[201221100] 


(1,0,2) 


[102112200] 


(2,0,2) 


[210112200] 


(1,1,1) 


[221120100] 


(2,1,1) 


[221101200] 


(1,1,2) 


[112102200] 


(2,1,2) 


[110212200] 


(1,2,1) 


[220121100] 


(2,2,1) 


[221201100] 


(1,2,2) 


[112202100] 


(2,2,2) 


[112210200] 



Table 2. 



4. De Bruijn Graph Homomorphisms 

In this section, we study homomorphisms between De Bruijn digraphs of differ- 
ent orders. In Subsection 14.11 we introduce and characterize homomorphisms between 
Bn+k{q) and Bn{q) for any integer k > 1 that perform like the homomorphisms of the 
previous section in the sense that taking the inverse by one such homomorphism of 
a De Bruijn cycle in -B„(g) dissects Bn+k{q) into adjacent vertex disjoint cycles. In 
Subsection 14.21 we single out the binary case which, due to its simplicity, admits a more 
concise characterization. 

4.1. Homomorphisms with Property (D). Recall that a graph homomorphism is 
a mapping that preserves the graph structure, so that an edge in Bn+k{q) is mapped 
to an edge in Bn{q)- The following proposition characterizes graph homomorphisms 
between De Bruijn digraphs of different orders. 

Proposition 4.1. A necessary and sufficient condition for a map Dn,k '■ Bn+k{q) 
Bn{q) to he a graph homomorphism is that 

Dn,k{^ = {dki^l-, . . . , a^fc+l), dk{x2, . . . , Xk+2), • • • , dk{Xn, . . . , Xn+k)) 

where x = (xi, . . . , Xn+k) and is any fixed function of k + 1 variables. 

Proof. Sufficiency is quite obvious so we will only prove the necessity. Let Dn,k{xi, . . . , 
Xn+k) be {xi, . . . ,Xn) where Xj = di{xi, . . . ,Xn+k) is a function from to Zg for 
alH = 1, . . . ,n. For all values of Xi, . . . , Xn+k, Xn+k+i {xi, . . . , Xn+k) is a predecessor 
of (x2, . . . ,Xn+k+i)- Hence, since Dn^k is a graph homomorphism the diagram below 
commutes, where the horizontal arrows indicate an edge in the De Bruijn digraph. 
That is, 

\Xl^ . . . , Xn-\-k) ^ \X2i • • • , Xn+ki ^^n+fe+lj 

Dn,k i [ Dn^k 

(•^1, • • • , Xn) > (-^2, • • • , -^n, -^n+l) 
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(4.2) 



di{Xi, . . . , Xn+k) — di-i{x2, • • • , Xn+k+l)] i — 2, ... ,71. 



To finish the proof we need to estabhsh that (1) di{xi, . . . , Xn+k) = dj{xi, . . . , Xn+k) 
for all i j and (2) di depends at most on x,, . . . , Xj+fc. We establish this by iterating 
Equations (14.21) for z = 2, . . . , n. To avoid confusion we will denote di by df and 
when it is applied to (xi, . . . , Xn+k) and (x2, . . . , Xn+k+i) respectively (the left and right 
sides of the diagram above). This is meant to remind us that, e.g., the first variable of 
(if is X2. 

First X2 = d2{xi, . . . , Xn+k) = df-{x2, . . . ,Xn+k+i), so that di does not depend on 
its {n + ky^ variable and d2 does not depend on its first variable. Next, x^ = 
d^{xi, . . . , Xn+k) = d2{x3, . . . , Xn+k+i) (uotiug that by the previous result df^ does not 
depend on its first variable X2). It follows that d2 does not depend on its (n + fc)*'^ 
variable and d^ does not depend on its first and second variables. 

Continuing this way we see that for i = 2, . . . , n, di-i does not depend on the 
last variable and di does not depend on its first {i — 1) variables. In particular, 
dn depends on at most x„, . . . , Xn+k and dn-i depends on at most Xn-i, • • • , Xn+k-i- 
Next dt _i{xn-i, • • • , Xn+k-i) — (-^fi-i ' • • • ' ^n+k) implies that dn-2 docs not depend 
on its (n + A; — l)*** variable. Continuing with Equations (14.21) iteratively and back- 
wards this time we establish requirement (2) above. But then Equations 14.21 reads 
df{xi, . . . , Xi+k) = df_^{xi, . . . , Xi+k). Hence, for all i di is a. fixed function of A; + 1 
variables which establishes (1). □ 

By Fact l2.1l a vertex disjoint cycle in 5„ is the D-morphic image of two vertex disjoint 
cycles in Bn+i starting respectively with zero and one. We generalize this D-morphism 
by finding and characterizing homomorphisms Dn^k so that each vertex disjoint cycle in 
Bn{q) is the image of exactly vertex disjoint paths in Bn+k{q), one for each starting 
string of size k. Such a homomorphism (and the corresponding function dk) will be 
said to have property (D). We will illustrate this property with some examples before 
we state Theorem 14.31 which characterizes functions dk that have property (D). In fact 
a direct inspection of the sixteen Boolean functions of two variables shows that the 
only homomorphisms with property (D) from Bn+i to 5„ are the D-morphism and its 
bitwise complement. 

Examples 4.2. (a) Consider the mapping D12 from B^ to Bi that uses the function 
d{xi, X2, X3) = Xi + X2. The inverse sets of and 1 are respectively {000, 001, 110, 111} 
and {010, Oil, 100, 101}. The edge (0, 1) of Bi is mapped back to the four edges 

{(001; 010), (001; Oil), (110; 100), (110; 101)} 

Note that, even though each edge in Bi is the image of four edges in B^, it is not 
possible to construct an edge starting with arbitrary strings of size two that is mapped 
to a given edge of Bi. For instance, there is no edge in B3 that starts with either 01 
or 10 and whose image is the edge (0; 1). 

(b) The function Hn^k from Bn+k{q) to Bn{q) for A; > and n > 1 was defined in 
[H] as Hn,k{xi, ■ ■ ■ , Xn+k) = {xk+i, " " " , Xn+k)- In othcr words, this function trims the 
k leftmost symbols of a word so as to make it a word of size n. Obviously, this is a 
homomorphism having, according to the notation of Proposition 14. 11 d{xi, ■ ■ ■ , Xk+i) = 

Xk+i, hence Theorem 14.31 below shows that it does not enjoy property (D). In fact the 
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inverses of any cycle in Bn{q) by H^^k disagree only in their first k terms while the 
body of the sequences are all equal to the original cycle. 

(c) Using the function X2, xs) = xi + X3 however, the edges (0; 0), (0; 1), (1; 0), 

(1; 1) of Bi are respectively mapped back to the following sets whose union constitutes 
the edge set of -B3, each edge appearing exactly once. 

{(000; 000), (010; 101), (101; 010), (111; 111)} 
{(000; 001), (010; 100), (101; Oil), (111; 110)} 
{(001; 010), (Oil; 111), (100; 000), (110; 101)} 
{(001; Oil), (Oil; 110), (100; 001), (110; 100)} 

Hence d*^^-* enjoys property (D) while d does not. 

Theorem 4.3. (a) A homomorphism from Bn+k{,<i) to Bn{q) that is induced by dk{xi, . . . , 
Xfc+i) enjoys property (D) if and only if dk is one to one in each of the variables xi and 
Xk+i when all the other variables are kept fixed. That is, if and only if dk{xi, . . . , Xk+i) 
defines a Latin square for each set of fixed values of X2, - ■ ■ , Xk- 

(b) The total number of homomorphisms with property (D) is {AqY where Aq is 
the number of q x q Latin squares. 

Proof. Part (b) follows directly from (a), see [TH] for more about the sequence Ag. To 
prove Part (a), first let dk be a function with property (D) and C = [ci,...,q] be 
an arbitrary vertex disjoint cycle in Bn{q). By definition of property (D) each word 
(xi, . . . , Xk) in can be appended by a symbol Xk+i so that dk{xi, . . . ,Xk, Xk+i) = Ci. 
This says that dk is surjective from to Zg (hence injective) with respect to the last 
variable. 

Now let x[ be such that dk{x'i, X2, ■ ■ ■ , Xk, Xk+i) and dk{xi,X2, . . . ,Xk, Xk+i) are equal 
to Ci. Since dk is bijective with respect to the last variable, there exist unique values 
Xk+2, • • • , Xn+k+1 such that Dn,k{x2, • • • , Xn+k+i) = (c2, • • • , c„+i). If x[ ^ Xi then the 
two distinct inverse edges 

, . . . , Xn+k] ^2, . . . , Xn+k+l)i (^^l; X2 ■ ■ ■ , X^+k] ^2, ■ ■ ■ ■, X^+k+l) 

share a common vertex, contradicting property (D). Hence dk is one-to-one in the first 
variable. This establishes the necessary condition. 

Conversely, let dk have the claimed form and let c = [ci, . . . , q] be a vertex disjoint 
cycle in -B„(g). Given any string xi, . . . , it is possible to find a value 6 G Zg so that 
dk{xi, . . . ,Xk,b) = Ci, since dk is surjective with respect to the last variable. Hence 
the value Ci has a set of q^ inverse images that includes all possible strings of size 
k as prefixes. The same argument implies that c has exactly inverse images. To 
show property (D) we need to show that no substring of size n + k occurs more than 
once in the collection of pre-images. Write the pre-images of c as a rectangular array 
(xjj); i = 1, . . . ,q'',j = l,...,k + l + n — 1 (corresponding to the "linearized" cycle 
ci, . . . , Q, ci, . . . , c„_i) where the set of prefixes of size k coincides with the q^ distinct 
words of this size. Let us denote by uJij{u) the substring of size u on the 2*'' row of 
(xij) that starts with Xij. Assume there exist integers ii, i2,ji,j2 such that ujij^j-^{n + k) 
coincides with oji^j^iji + k). Obviously ji ^ j2 implies that a string of size n occurs 
twice in c, thus contradicting the assumption that it is vertex disjoint. Assume then 
that jo is the smallest integer with u;j^j,)(n + k) = uJi,^^jg{n + k), which in particular 
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means that Ui-^jg^k) = LJi^j^^{k). By construction of (xij), jo > 1- Hence ^^(tUij + 
1)) = '^fc('^j2,jo-i(^ + 1)) = Since the last k components of cUj^ + 1) and 

C(jj2j„_i(A; + 1) are the same, the one-to-one-ness of dk with respect to the first variable 
implies that = Xi^j^-i. Therefore ci;j^jQ_i(n + k) = Ui^jg-i^n + k), which 

contradicts the minimality of jo- This establishes the theorem. □ 

In general we see that applying the inverse of a homomorphism to an vertex disjoint 
cycle in Bn{q) creates multiple cycles in Bn+k{.<i)- If Bn{q) is partitioned into vertex 
disjoint cycles then the inverse homomorphism naturally induces a partition of Bn+k{(l) 
into vertex disjoint cycles. 

4.2. The Binary case. We treat here the binary case separately because its simplicity 
allows for a more concise characterization of the shape of homomorphisms with property 
(D). 

Theorem 4.4. A necessary and sufficient condition for a homomorphism Dn,k from 
Bn+k to Bn to have property (D) is that 

dk{xi, Xk+i) = xi + h{x2, ...,Xk)+ Xk+i, 
where h{x2, ■ ■ ■ , Xk) is any Boolean function of k — 1 variables. 

Proof. By Theorem 14.31 we only need to show that a binary function dk is bijective 
with respect to the first and last variables if and only if it has the form claimed in this 
Theorem. In effect, if dk{xi, . . . , Xk+i) is bijective in xi and in Xk+i then it satisfies the 
equations 

dk{xi,X2, . . . , Xk+l) = 1 - dk{xi,X2, . . . , Xk+l) = dk{xi,X2, . . . , Xk+l). 

So that dk{xi, X2, . . . , Xk, Xk+i) = dk{xi, X2, . . . , Xk, Xk+i). 

Therefore for each fixed set of values for X2, ■ ■ ■ , Xk, dk{xi, . . . , Xk+i) = dx2,...,Xk{^i,Xk+i) 
is either xi + Xk+i or xi + Xk+i + 1. This can be rephrased to establish the necessity. 
The converse is obvious because dk is linear in the first and last variables. □ 

This elegant form of dk is mainly due to the "lack" of terms in Z2. While Theorem l4.3l 
shows that dk = OiXi + h{x2, ■ ■ ■ , Xk) + [3xk+i is sufficient for property (D), the following 
example illustrates why property (D) homomorphisms can not be all written in such a 
simple form even for g = 3. In fact all the twelve 3x3 Latin squares can be written 
in function form as 62) = tti&i + 02^2 + 03 where 6,, G Z3, ai 7^ and ^2 7^ 0. 
Fr the values 0, 1, 2 of X2 let dk{xi,X2, xs) be respectively xi + X2 + X2, 2xi + X3 and 
Xi + 2a;2 + 2x3. Then dk has property (D) by Theorem 14.31 but it is not linear in either 
Xi or Xk+l, despite the simple form of Latin squares. Notice that when g > 3 most 
Latin squares are already nonlinear. 

While the only binary homomorphism for = 1 is Lempel's D-morphism (and its 
bitwise complement), there are essentially two homomorphisms for k = 2 that are 
induced by the functions d^^^ = X1+X3 and d^'^^ = X1 + X2 + X3. Note that the former is 
just the D-morphism iterated twice. The only other two homomorphisms are bitwise 
complements of d^^^ and d^'^\ The cases k > 3 allow for nonlinear homomorphisms 
such as d{xi, . . . , x^) = Xi + X2X3 + X4. 

Let c = [ci . . . Ci] be an arbitrary but fixed cycle in Bn, started at a fixed word, say, 
. . . 0. Then for each homomorphism Dn,k with property (D), c defines a map Vc on 
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the set Z2 as follows. T>c{zi . . . Zk) is the suffix of length k of T)^ started at the string 
Zi - . . zi^. The inverse image is generated by the recursive equation 

Zi = d-k + Zi^k + h{zi-k+i, . . ., Zi-i); i = k + 1, . . . ,k + I, 

where h is as in Theorem 14.41 and zi, . . . ,Zk are the required initial conditions. It can 
be seen that property (D) implies that Vc is a bijection. When the D-morphism is 
used, any De Bruijn cycle b„ yields the identity permutation on the set {0, 1}. This is 
a restatement of the fact that the inverse image of any De Bruijn cycle b„ under the 
D-morphism makes two dual cycles in Bn+i- Since a binary De Bruijn cycle necessarily 
has an even number of ones, this follows immediately from Fact 1 above. The next 
proposition concerns the function S'^^ defined above. 

Proposition 4.5. For any integer n > 1 and any De Bruijn cycle b„ = [61 . . . 621], the 
homomorphism induced by the Boolean function d^'^\xi, X2, X3) = xi+ X2 + X3 defines a 
permutation of the set of seeds {00, 01, 10, 11} with exactly one fixed point Z1Z2 obtained 
by 

Zi = aQ + Sfi^ai + 5n,ici2j Z2 = cii + Snfia2 + 6n^iaQ, 
where n = n mod 2, aj = a^"'' := Yli^3i+j i^od 2;j = 0, 1,2, the sum is taken over 
the range of indices of b„ (I < Si + j < 2"-), and addition in the index of aj is taken 
modulo 3. 

In other words, exactly one of the four sequences that form the preimage of b„ is a 
closed cycle in Bn+2- As a result, the other sequences together form one cycle of length 
3 ■ 2". 

Example 4.6. Let = [00011101]. We see that n = 1, Oq = 1, ai = 1 and 02 = so 
that the fixed point is Z1Z2 = 10. Indeed the inverse image by d^'^^ gives the following 
four sequences. 

0000010001; 0110100111; 1011001010; llOlllllOO- 

So the fixed point gives the only cycle of length 8 while the other three sequences make 
the following cycle of length 24: 

[oooooioooiioioomoiiiii]. 

Proof, (of Proposition 14. 5p Let i = i mod 3. Iterating the relation Zi = bi-2 + Zi^i + 
Zi_2, which is satisfied by the sequence {-Zij^^s ^, we get (for all i in the range of the 
latter sequence) 

LV3J-1 

^ (bsj+i + bsj+i+i) + ki^i + ^l2Z2 + Sio{zi + Z2). 

j=0 

where we define 60 to be zero. Note that 2" mod 3 = 1 or 2 when n is even or odd 
respectively. For each of these two cases, using the above recursive equation and the 
requirement Z2n+j = zj; j = 1,2 yields two linearly independent equations whose 
unique solution is as claimed. □ 

Shifting hn by a number that is not a multiple of 3 permutes the numbers aj, j = 
0,1,2. So it changes the permutation but still keeps one fixed point. This result is 
interesting because it is independent of the De Bruijn cycle used. The permutation 
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induced by may or may not have a fixed point, depending on b„. As a result, 
The function S"^^ can be used to generate De Bruijn cycles recursively by joining the 
shorter cycle (the one started at the fixed point) to the long cycle made of the other 
three starting values. Proposition 14.51 describes the way to identify the two starting 
digits of the short cycle. There is no simple way to identify a pair of conjugate words 
to perform this cross-join operation a priori, for example the alternating strings may or 
may not be on the same cycle. Note that the existence of a word without a conjugate 
on the shorter cycle is guaranteed because otherwise the cycle must be a De Bruijn 
cycle, see [9j. Consequently, one can find a cross-join pair by only searching the shorter 
cycle for a word without a conjugate there. This search takes 0{N'^) in the worst case, 
where = 2" is the length of the short cycle. This is manageable for small to medium 
word size n. 
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